% this codes the table for the function in Kim et al. (2010) to be
% used in the computaton of reservation values

% clear all
% clc

fun = @(zeta) normpdf(zeta) - zeta.*(1-normcdf(zeta));

% NB: input of fun  = r/sigma_z
%     output of fun = c/(beta*sigma_z)

input_list = linspace(-15,15,100000);

output_list = fun(input_list);

% plot(zeta_val,x_val)





